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The development of a speed-of-light hard-kill weapon system for military 
applications represents a significant advancement in technology over present 
conventional kinetic weapon systems. Over the past two decades, the US Navy has 
successfully developed a megawatt-class chemical laser; however, under some maritime 
environments, the high power beam propagation was unable to delivery sufficient energy 
to kill a modem anti-ship missile (ASM) due to significant atmospheric absorption and 
the resulting thermal blooming process. A critical problem to resolve for the shipboard 
high-energy laser weapon systems is to develop a shipboard-compatible megawatt-class 
laser weapon at a wavelength where the atmospheric absorption is smallest. The 
megawatt-class Free Electron Laser (FEL) has significant advantages over conventional 
weapon systems and other chemical high-energy laser systems. Infinite magazine, rapid 
response, and wavelength tunability make the FEL a suitable and desirable shipboard 
weapon system. 

This thesis divides into four chapters. Chapters I and II introduce the FEL and 
background theory of the FEL. Chapter III explores the analysis of the LANL 
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1. INTRODUCTION 



A. DEFENSE AT THE SPEED OF LIGHT 

In today's technology, weapon systems rely predominantly on kinetic kill, striking 
a missile with another missile or with several bullets. These weapons are relatively slow 
and require a longer detect-to-engage sequence when factoring in tracking time, launch 
time and time of flight to intercept. In searching for faster and faster weapons and ever- 
decreasing response times. Directed Energy Weapons (DEW) have been on the U.S. 
military's "love to have" list for several decades. DEWs are attractive for several reasons, 
and the most predominant reason is that they offer speed-of-light defense. 

In 1971, the U.S. Navy established the Navy High Energy Laser (HEL) Program 
with the intentions of developing speed-of-light hard-kill weapon system. Initial 
development focused on building a CO 2 gas dynamic laser technology, but a few years 
later its emphasis shifted to continuous-wave deuterium fluoride (DF) chemical lasers 
because they propagate far better in the maritime environments. In the early 1980s, the 
megawatt class Mid-Infrared Advanced Chemical Laser (MIRACL) was built by TRW 
for the Navy as the highest average power continuous wave (CW) laser in the west. The 
SeaLite Beam Director (SLBD) was built by Hughes Aircraft to accommodate and direct 
the MIRACL beam. The MIRACL and SLBD have successfully tested at White Sands 
Missile Range, New Mexico, in engaging and destroying several types of missiles in 
flight [1]. 

Although the MIRACLE/SLBD system was successfully tested against crossing 
targets in a fleet defense scenario, its lethality was not sufficient in some ship self-defense 
engagement scenarios due to insufficient crosswind [1]. Insufficient crosswind will cause 
the laser beam heating the air in its own path resulting in a defocusing optical lens known 
as thermal blooming [2]. The thermal blooming effect can severely reduce the High 
Energy Laser Weapon Systems' (HELWS) lethal range. As a result, it's necessary to 
replace the MIRACL with a megawatt-class laser device at a wavelength where 
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atmospheric absorption is smallest. An ideal solution is the megawatt-class Free Electron 
Laser (FEL) that is tunable to the optimum propagation wavelength. 

B. SHORT HISTORY OF FREE ELECTRON LASER 

The first generation of coherent radiation from free electron beams was 
microwave tubes that relied on slow-wave structures, and their development received a 
significant impact from the radar development during World War II. In 1951, Motz 
showed that an electron beam propagating through an undulator magnet c£in be used to 
amplify radiation. After 1 960, research on short wavelength lasers has been dominant, but 
they require an atomic or molecular medium to operate, resulting in a limitation on 
wavelength tunability. In 1971, J. M. J. Madey of Stanford University brought attention to 
what he proposed as the 'free electron laser’. Relying on research of synchrotron radiation 
sources, Madey conceived a device that would lase in the visible using a beam of 
relativistic electrons as the source. In 1976, the first successful free electron laser (FEL) 
experiment was demonstrated by Madey and his co-workers at Stanford University, in 
that they measured 7% gain from an FEL configured as an amplifier at 10pm 
wavelength. One year later, another successful FEL operation of the same FEL configured 
as an oscillator. These successful experiments created a large interest in FEL research. 
Since then, scientist and experimental groups around the world built FELs that operate at 
frequencies ranging from microwave to the UV [11]. Three other FELs were built soon 
after. The first was built at Laboratoire pour I'Utilisation du Rayonment 
Electromagnetique in Orsay, France, which operated in the visible range. The second 
involved a team from TRW at Stanford using the Superconducting Accelerator (SCA), 
lasing in the near infrared (IR). The third was built at Los Alamos that lased in the mid- 
infrared [3]. Now, there are many FELs used for research around the world. 
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c. 



BASIC FREE ELECTRON LASER THEORY 



The FEL is a device that transforms the kinetic energy of a relativistic electron 
beam into electromagnetic (EM) radiation [4]. This is accomplished by sending a 
relativistic electron beam from a linear accelerator, microtron or a storage ring etc., 
through an alternating magnetic field produced by a device known as an undulator as 
shown in Figure 1.1. The interaction of the negatively charged electrons with the spatially 
periodic magnetic field induces a transverse acceleration on the electrons that causes 
them to exchange energy with a copropagating radiation field. The kinetic energy 
extracted from the electron beam is converted to electromagnetic radiation. The energy of 
the electron beam and the period of the undulator magnet determine the wavelength of 
radiation. This system provides a powerful source of tunable, coherent, electromagnetic 
radiation that can operate in a wavelength range from about 1 cm to the X-ray region. 
This system can also provide continuous wave (CW) operation, or very short pulses in the 
nanosecond, picosecond, or even subpicosecond region. 



Electron Beam 




Figure 1.1. Schematic diagram of a free electron laser. 



3 



4 



II. PRINCIPLES OF OPERATION OF FEL 



A. 



ELECTRON TRAJECTORY IN THE UNDULATOR 



For better understanding of the physical mechanism underlying the operation of an 
FEL, we first consider the interaction of a free electron in the undulator and an 
electromagnetic wave. The Lorentz force equation and the equation for the energy 
exchange between the electron and the electromagnetic wave are 



P # I 

F = —=-^ 
dt 



V 

E+ — X 5 
c 



d(ymc^) 



=-ev-E, 



( 2 . 1 ) 



( 2 . 2 ) 



where e, m and v are the charge, mass and velocity of the electron, respectively; y is the 
relativistic factor of the electron, y= l/-^(l ; y0=v/c; p is the relativistic 

momentum of the electron, p=ymv ; E is the electric field of the EM wave; B is the 
magnetic field ; ymc^ is the electron energy; and c is the speed of light in vacuum [5]. 

Considering an electron with velocity v = 0 c and a transverse electromagnetic 
wave propagating colinearly, in the z-direction in a vacuum, then equations (2.1) and 
(2.2) become 

d{y 0) e / = 

dt me 



;(£ + /?x5). 



(2.3) 



dy e _ _ 
-j-= — p-E. 
dt me 



(2.4) 



Assume the undulator has a helical magnetic field of the form 5 = 5[cos(A:^2),sin(^p2),0], 
where k^=2n ! is the undulator wavenumber, and \ is the undulator wavelength. 
The corresponding circular-polarized, plane-wave optical field is then 
= E (cos sin ^,0), and B^=E(sinip,eosy/,0) , 
where y/=kz-cot + ^ . £ is the electric and magnetic field amplitude in cgs units, t is 
time, ^ is the optical phase, k = col e is the optical wavenumber, and O) is the optical 
angular frequency. Substituting the fields into Lorentz force equations yields 
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d{yP\) ^ 

— T — =- — [£(l-y9,)(cosv/,-smv/,0) + /?, B(-sin(^^z),cos(A:„2),0)], (2.5) 
at me 

d{yP,) e 

— ^ = - — [E{P^cosi{/- p + B{p^sm{k^z)- p cos.{k^z))\, (2.6) 

at me 

and ^=--^E{p ^cosi!/- PySini!/), (2.7) 

where P = Pi_ + P., P ^ = P p + Pyj , and p. = p,k. For relativistic electrons, /?, «1, 
hence £ (1 - p,) « p, 5, so in the right-hand side of (2.5) the first term is ignored. The 
second term is 

sin( k^z),cos( k^z),0) . (2.8) 

dt me 

Integrating the transverse equation (2.8) and assuming a perfect injection of electron 
beam into the undulator cavity, which makes the constant of integration zero, yields 

rPL=-{-^^^(cos(k^z),sin(k^z),0 ) . (2.9) 

Divide both side of (2.9) by y to get 

P^i^)= (cos(A:„z),sin(^„z),0), (2.10) 

7 

where the undulator parameter is K=eB / k^me^ =eBA^ / 2mne^ . Equation (2.10) 
describes the transverse motion of electrons in the helical undulator [6]. 

B. FEL RESONANCE - ELECTRON AND PHOTON RACE 

Accelerated electrons traveling at a relativistic speed along the longitudinal axis of 
the undulator of the FEL radiate optical energy [6]. In order for the FEL to have gain (i.e. 
the fractional change in optical power in a single pass through the undulator), a net 
transfer of energy from the electron beam to the optical beam must occur. In order for an 
optimum exchange of energy from an electron to the optical field to occur, the electron 
must oscillate in phase with the optical field. This optimum energy exchange occurs when 
the electron falls behind the optical fields by a distance of one optical wavelength within 
the span of one undulator period [6]. 
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speed of electron P 2 ^ 




Speed of light c 




Light wavelength 



t=0 Undulator wavelength t=4,/ p 2 ^ 

\ \ 







Figure 2.1. Electron-Photon race. 



Consider an optical wave that travels at speed c and passes over an electron traveling at 
speed Pf as shown in Figure 2.1. At resonance condition, one wavelength of the optical 
wave passes over an electron while it travels through one undulator wavelength so 



that light wavelength 2,= (c - P^c)t . Insert t= X^l P^c yielding X=- 



\_zA 



Xg. Substitute 



Pi = — into = \- P^ = Pi - Pi yields 

y 

{\ + K^)y-^ = \-Pl . ( 2 . 11 ) 
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\ + K^ 



Since » 1 , 1 - P. 



2r‘ 



, then 



/l = 



1 + /: 
2y^ ' 



( 2 . 12 ) 



The equation (2.12) is the resonance condition that demonstrates the tunability of the 
PEL. By adjusting the undulator parameter K, or the energy of the electron beam 
(y-\)mc^ or the undulator wavelength the wavelength of the optical field can be 



tuned to a desired wavelength within a prescribed bandwidth [7]. This is a significant 
advantage over conventional lasers that are typically constrained to a fixed wavelength. 



C. PENDULUM EQUATION 



This section derives the equation of motion of the electrons resulting from the 
interaction with the optical field. Substitute (2.10) into (2.7) to get 

eKF 

r = —cos(C + ^), (2.13) 

ymc 

where 

^ = I// + k^z = {k + kg)z - oX + ^ . (2.14) 



For -n !2< ^ + <f> <n 12, cosine function in (2.13) is positive and the electrons will 
gain energy from the optical field, whereas for 7 t ! 2 < <^ + (j> ! 2, the cosine is 

negative and the electrons will lose energy to the optical field. Differentiating (2.1 1) with 
respect to time t gives 2y~^y(l + )= 2ft, p, . Relating y to p, gives 



f _ y^PzPz 

y (1 + ^^)’ 



(2.15) 



Differentiating (2.14) once yields ^ = (k + k^)z-o) = (k + k^)cp, -co, and 
differentiating twice yields ^ = (k + k^)cp,. Insert the resulting expression for/?^ into 



(2.15) to get 



y ^ y^Pz^ 

y {\ + K^){k + ky 



(2.16) 
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Since X « then « k, so {k + k^)c^kc = a>. Using this approximation and the 
resonance condition (2.12), equation (2.16) becomes 

y (} + K^)o} IcOg 



So, 




y eKE 
y y^mc 



cos((^ + <!>), then the electron equation of motion is 



^■ = — \ cos{C + (l>). (2.18) 

y me 

Equation (2.18) describes the electron phase dynamics in the form of the pendulum 
equation [6]. 

Consider time evolution through the undulator L / for an average electron, 



where L is the length of the undulator. Introducing relevant dimensionless time 
T = PfjCt I L » ct I L , v/e see that r will vary between 0 and 1 along the undulator length. 
Substituting dt = (Z, / c)r/r into the pendulum equation (2.18) gives 



C = v = \a\cos{C + (!>), (2.19) 

o 

where the notation {..) = d{..)l dr , \a\ = AnNeKLE I y^mc^ is the dimensionless optical 

field strength, and N = L! is the total number of undulator periods. Equation (2.19) is 

the electron pendulum equation that governs the phase-space motion of the electrons 
within the undulator and under the influence of the optical wave [6]. 
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D. 



OPTICAL WAVE EQUATION 



An electromagnetic wave will propagate in the absence of a current density 
according to 

\ 



2 1 



V" - 



c' de 



A{r,t) = 0, 



( 2 . 20 ) 



where the Laplacian is + 5^ , r is the vector position, t is time, and^i is the 

optical vector potential [5]. The optical electric and magnetic fields can be expressed in 
terms of the vector potential by the following relationships 



= 1 dA 

E = — , 

c dt 



( 2 . 21 ) 



and B = V^A. (2.22) 

Assuming that the optical beam is a circularly polarized plane wave, then the vector 
potential has the form 

A = — ^^^(sinvji,cosvj/,0)=9?c , (2.23) 

where S ' , y/ = kz - aX + <j> is the plane-wave phase, e = (-/,1,0) is the 

polarization vector, k = 2nd (o is the optical wave number, a> is the optical frequency, 
and (!>{r,t) is the optical phase [7]. For convenience, we write 



A = jS{r,t)d'“ , 
k 

where a = kz-cot . Differentiating (2.24) with respect to t gives 



(2.24) 



dA e(^ . 
— = —\® -1 
dt k^ 



\co 



(2.25) 



where /dt . Differentiating (2.25) with respect to t gives 



d^A 



= — {s - 2x0) S - 

dt^ k^ 



O) 



(2.26) 



where / dt^ . Differentiate (2.24) with respect to z gives 
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(2.27) 



dz ^ 

where (..)'=d(..)/ dz . Differentiate (2.27) with respect to z gives 

= |(J’" + 2ik , (2.28) 

where =d^(..)/ dz' . Since the laser beam consists of a narrow range of frequencies, 
the variation in the electric field amplitude and phase with respect to space over an 
optical wavelength A is small, so that k S’ and <{>' The variation in the electric 

field amplitude and phase with respect to time over an optical period is also small, so that 
«co ^ and [7]. Therefore, (f"is negligible compare to k^' , and ^ is 

negligible compare to co ^ , then (2.26) can be simplified to 

dt^ k ^ ’ 

and (2.28) can be simplified as 

^ = -^(2ik^'-A:'(r). 

Substitute (2.29) and (2.30) into the left hand side of (2.20) yields 



1 a' ^ 


A(r,t)= 


V\^ + lik 




c dt ) 


k 




1 ^ J. 



(2.29) 



(2.30) 



where V , = — ^--i- 






. Hence, 







1 d\ 

V 


V' +2ik 


H 


1 

1- 


i,,5z 


c dt ) 



= 0 . 



(2.32) 



The first partial second derivatives with respect to jc and y describe how the beam changes 
due to diffraction as it propagates. The last partial derivatives with respect to z and t 
describe the propagation of a plane wave along the z axis. 

Introduce dimensionless time r-ctIL so that r varies from 0 to 1 along the 
undulator length. It's convenient to follow the wavefront with a new coordinate 
'z = z-ct . With this coordinate change. 



11 



(2.33) 



dz 



(dz\ 



dz 



\U^ J 



5z 



• + 



A 

^dzjdr dz' 



since dz / dz = \ and dz/ dz - 0 . Also, 



lA-i 

c dt c 



^ dz^ 
\dt j 



dz 



ydt ^ 



dr 



__A lA 

dz^ Ldz 



(2.34) 



since dz/dt = -c and dz/dt = c/ L. Therefore, the propagation operator in (2.32) 
becomes 



A lA-A 

dz c dt dz 



r A 1 ^ ^ Id 

dz L dz . 



J 



L dz 



(2.35) 



Hence, the full wave operator is 



V' +2/1 







V^+2/1-— 


_L 


^dz c dt J 







(2.36) 



Introducing the dimensionless transverse variables x = x^fkTlL and y = y^JkilL , so 

2 2 2 2 

that -^-Y = {k/2L)^-r^ and -^ = (A:/2L)-^ . Hence V\={k/2L 
dx dc' dy ^ 






(A:/2I)V^, where = 



dc^ 



. Substituting into (2.36) and (2.32) yields 

= 0. (2.37) 






4 dz_ 

Equation (2.36) is known as the parabolic wave equation [7]. In dimensionless 






where 



parameters, the optical field is represented by a = \a\e’ 

|a| = 4;r NeKLE I {y^ mc^ ) cc E , so W and a differ by a multiplicative constant. Therefore, 
exchange ^ for a, the dimensionless parabolic wave equation (2.37) becomes 






a = 0. 



(2.38) 



4 dr. 

In the presence of the electron beam current density, the optical wave equation 
(2.20) is 
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2 1 5 



2 \ 



V" - 



c' dt^ 



Ak - 
A{r,t) = J^, 



(2.39) 



where J ^ is the transverse electron beam current density [8], Substituting (2.30) into 
(2.39) yields 

ee'‘ 



2 ,( ^ 


1 


V\+2ik\ — 






c dt J_ 



_ An -i 
c 



(2.40) 



Performing a dot product with e =(/,l,0) on both sides of (2.40) and using 



e’ - e = (U,0) • (-/,1,0) = 2 , yields 

^ e~'“ . (2.41) 

c 

The total transverse electron beam current is the sum of all single particle currents. The 
single particle current for one electron in the beam is = -ecp^5^^\x -r^) [8] where 

is the position of the /'* electron and -^) is the three-dimensional Dirac delta- 

function. The total beam current is 

J, (2.42) 

I t 

From (2.10), the transverse motion is p^={- K / y\coskgZ,sink^z,{))= 
l(-/i:/y)/ec For convenience, we write = — —iee , so that 

r 

-i(kj,z+a) 

-eV'" = iecK{e' •e)T ^ S^^\x-r ,) . (2.43) 

. y 

We can simplify the sum ^ by assuming a local electron density p, so that (2.41) 

I 

becomes 



1 9 

-V]+ik 

2 ^ 



1 

\dz c dt 



-V{ + ik 
2 ^ 



dz c dt 



^ = -2nieKp 




(2.44) 
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where 4" = + a is the electron phase in the combined undulator and optical fields, 

and (...) is a weighted average of sampled electrons at a specific site which replaces the 
summation over all electrons within one optical wavelength. Define the dimensionless 
beam current j = %N{e7r KLY pl{y^mc^) , the wave equation becomes 

a(x,y,z,z)=- j{e~‘^Y (2.45) 

In some cases where the optical wave diffraction is not significant, we can ignore V ^ in 
(2.45), then 




°a = -j{e-'^). (2.46) 

The optical wave equation (2.46) couples the electrons to the copropagating optical wave 
by relating the electron phase to the change in the optical field. The equation (2.46) can 
be expressed in terms of amplitude and phase. 



H = -y(cos(4' +</>)), 

and ^ = |^(sin(4' + ^)) • 



(2.47) 

(2.48) 



Equation (2.47) reveals that when the electrons become bunched near + <p)^ ^ , the 
optical field amplitude will grow, and when bunched near (4” + ^) « 0 , the optical field 
amplitude will decrease. Equation (2.48) describes the optical phase evolution. 
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E. GAIN DEGRADATION DUE TO ELECTRON BEAM QUALITY 



In the design of an PEL, there is often a trade-off between high beam current and 
good beam quality as determined by the accelerator. The dimensionless current 
j oc where 1 is the actual beam current , N is the number of undulator periods, 

and X is the optical wavelength. Usually the optical wavelength X is fixed for a given 
application so that j is maximized by adjusting current I and the undulator periods N. 
Increasing / tends to degrade the beam quality from the accelerator, whereas increasing N 
tends to increase the PEL sensitivity to the beam quality because of narrowing the gain 
spectrum bandwidth [9]. 

Emittance e describes one aspect of electron beam quality. A real electron beam 
has an emittance e describing how electrons enter the undulator with a range of angles 
9^ with respect to the undulator axis, and at a range of transverse positions slightly 
off axis. The electron beam may also have a small energy spread txyly . The energy 
spread and emittance of the electron beam translate into a spread in phase velocities in 
phase-space. The phase velocity spread due to energy spread is given by [10] 

ts.v =AiiN t^y ! y , (2.49) 

and the phase velocity spread of an electron beam due to emittance is given by 

1 H- A 

Electrons enter the undulator with a spread in phase velocities due to a spread in 
energies, and are randomly distributed in both the longitudinal and transverse directions. 
The typical distribution of electron phase velocities is assumed to have a Gaussian form 
about initial phase velocity , 



/o(?) = 









(2.51) 



where <Jq = An N ls.y ! y is the standard deviation and q is the phase velocity away from 
distributed [9]. The characteristic function P’o(r) of the distribution /c(^)is defined 



as 
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Evaluating the above integral equation for a Gaussian yields 

If the beam quality is perfect so that cTq = 0 , then Fg (r) =1 . However, when the beam is 
not perfect, the characteristic function |E’(j(t)| decays over the length of the undulator in 
a characteristic time on the order of \! [9]. 

Poor beam quality reduces the single-pass FEE gain by degrading electron 
bunching. As an example, consider an FEE design with N = 25 periods having a 
dimensionless current density of / =10 and initial phase velocity of =3 operating with 

an initial dimensionless optical field a<, = 1. In Figure 2.2, electron beam quality is 
assumed to be nearly perfect, i.e. = \ . Some electrons gain energy and move ahead of 

the average flow, while other electrons lose energy to the radiation field and move back 
behind the average flow causing visible spatial bunching. 



r *** PEL Phase Space Evolution ***’“ 
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Figure 2.2. Phase space evolution for a nearly perfect electron beam quality. 
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In Figure 2.3, the longitudinal energy spread increases to <7^=6, and this 

degradation in beam quality causes much less bunching when comparing to Figure 2.2. 
The poor beam quality is responsible for a decrease in final gain to G/= 0.10534 and a 
smaller optical phase shift <l> . 
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Figure 2.3. Phase space evolution for low gain with electron beam degradation. 

F. TAPERED UNDULATOR 

The characteristics of the FEL interaction can be altered by changing the 
undulator properties, i.e. undulator polarization, wavelength, or field strength, along its 
length [5]. As an FEL reaches saturation, the electrons have given up enough energy so 
that they are no longer in resonance with the optical wave, which decreases the gain. 
Tapering changes the imdulator's resonance condition along its length so that the 
electrons remain in resonance longer, hence increasing the overall gain in strong fields 
[7]. From the resonance condition (2.12), in order to extract energy out of the electron 
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beam, i.e. reducing y, and convert into radiation field and still maintain resonance 
condition, we can either taper the undulator wavelength , i.e. reducing A^, or taper the 

undulator parameter K, i.e. reducing the undulator field strength B. However, tapering K 
is easier since the undulator gap needs to be increased to decrease B. Tapering essentially 
has the same mathematical result as accelerating the electrons. This artificial acceleration 
is represented by a dimensionless parameter S. A tapered undulator represents an 
advantage over the untapered undulator design when the artificial acceleration S exceeds 
the deceleration that can be obtained without taper in strong optical fields. When the 
undulator wavelength is decreased, the artificial acceleration S is given by 
S w -2kNAX^ / A„ . When the undulator field strength is decreased, 

5 fs-47iNK^ABIB{\ + K^). With tapering is included in the undulator, the pendulum 
and wave equations become 

I = v = 5 + \a\cos{i; + <f>), ( 2 . 52 ) 

O 

and a = -j < e~'^ >, 

where ^ = ^k^{r')dr' + kz - cot [10]. The wave equation remains the same, but the 

simple pendulum equation includes an additional constant torque due to the artificial 
acceleration S . 

I I I 1^/2 

The necessary and desirable criteria for tapering is |a|><5>4p| > Itt . The 

condition on the left is for trapping electrons in the tapered phase space; the middle 
condition is that the tapered acceleration exceeds the natural untapered deceleration; and 
the condition on the right is that the tapered PEL works in the strong field regime [10]. 

Electrons initially near the phase ~0 are accelerated by both the strong optical 

field |a| and the artificial "torque". The electrons initially near «;r undergo a 

cancellation between the optical field and the torque, leaving them trapped in closed 
orbits. The taper is effective because electrons near the phase for gain are trapped, 

and continue to contribute to the interaction. Electrons near the phase for absorption 
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4" »0 are taken away from resonance and eventually stop interacting. The imbalance 
leads to a net gain [10]. 

Compared to the untapered undulator, the tapered undulator is more efficient in 
strong optical fields, but has smaller gain in weak optical fields. The efficiency of a 
tapered undulator is estimated at rjg =S l{S;r N) . For example, if ^»280;r for an 
undulator with 23% taper in wavelength over N = 300 periods, then the efficiency is 
estimated at rjg f>ilQ/(2N ) which increases by a factor of 70 over the natural efficiency 

rj' fs\t(2N) [10]. Tapering is a design to obtain higher overall gain and efficiency in 
strong optical fields. 

An example of a tapered undulator in strong optical fields is shown in Figure 2.4 
with an untapered undulator shown in Figure 2.5. Both designs are based on an undulator 
length of I = 6 m corresponding Xo N= 300 periods. The dimensionless current density is 

j = 12000 y>n in the high gain regime, and initial field =150 »7t in the strong optical 

field regime. In Figure 2.4, the tapering rate S = 260;r turns on at = 0.19 along the 
undulator resulting in a final gain of Gf = 556 and FEL extraction efficiency of 
ri«12.5%. Initially the gain increases exponentially until =0.19 where the taper 

turns on and saturation is reached. With tapering, the gain continues to increase slowly 
instead of decreasing after saturation, as shown in the upper right-hand portion of the 
Figure 2.5. The phase velocity of the untrapped electrons near v «740 is shifted by the 
tapering acceleration S whereas the trapped electrons remain near resonance due to 
strong optical fields. The lower right-hand portion of the Figure 2.4 is the optical phase 
shift, and this phase shift A(/> « 1 .4;r over the length of the undulator is a characteristic of 
high gain. 
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*** FU Phase Space Evolntion *** 
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Figure 2.4. Phase space evolution for high gain in strong fields with a tapered undulator. 
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Figure 2.5. Phase space evolution for an untapered undulator. 
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III. SIMULATIONS OF THE REGENERATIVE 
MW FEL AMPLIFIER 

A. MW-CLASS FEL CALCULATIONS 

In order to destroy an incoming missile, the laser power required on the missile 
surface is about 10 kW/cm over a 5 cm radius spot for a two second duration. The 

extinction coefficient due to aerosols at sea level is a « 0.03 km’’, and through e~‘^ 
describes the absorption loss at X « 1 p m wavelength over a distance z. Let P, be the 
Transmitted Power that must leave the ship to destroy a missile at a range P, then the 
power at the target is -P, . If the range of the missile is 5 km, the laser power 

which must leave the ship is then P, = [10 kW/cm^] [7i(5 cm)^] « 910 kW. 

As an example, consider an FEL with N - 300 periods using a beam of 100 MeV 
electrons with a peak current of 400 A, then the peak electron beam power P = 400 
Ax 100 MV « 40 GW. The micropulses generated by the FEL have a length of 0.001 m 
and are separated by 6 m, then the fraction D = 0.001/6 denotes the duty cycle of the FEL. 

The average power generated within the FEL is given by P=PxD»6.7 MW. The 
efficiency necessary to destroy an ASM is r| = P, /P»^13%. Therefore, the single-pass 
FEL efficiency r| needed to supply the Laser Power to destroy the incoming missile is 
13%. An approximation of the natural single-pass efficiency of an FEL is only 
T] = 1 / 2N = 2% so that a tapered undulator with high efficiency is needed. The tapered 
undulator described in the previous section would be more than adequate. This type of 
illustration is the essential driving factor in determining the FEL parameters for possible 
MW-class weapon systems. 
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B. LANL REGENERATIVE MW EEL AMPLIFIER 



The US Navy is investigating the possibility of using a MW-class Free Electron 
Laser for ship self-defense against anti-ship missiles (ASM). A design workshop resulted 
in two possible MW FELs, the oscillator and the regenerative amplifier. The MW 
oscillator design, although more compact overall, presents a challenge in recirculating an 
intense electron beam because of a phenomena known as coherent synchrotron radiation 
(CSR) feedback [11]. The regenerative amplifier design uses a longer undulator and 
relies on higher extraction efficiency to achieve high average power so that less current is 
needed from the accelerator and no recirculation is required. The regenerative amplifier 
design has been proposed by LANL and Boeing with the analysis being performed at 
NPS. This section explores the analysis of regenerative amplifier design and optimizes its 
efficiency. 

A simple schematic of the regenerative amplifier design is shown in Figure 3.1. 




Photo Injector 



Feedback Loop 




Electron Dun^ 



Figure 3.1. Conceptual MW class Amplifier FEL. 



The amplifier design does not use energy recovery so there is no requirement to 
bend the high energy beam which can cause deterioration of the electron beam quality. 
However, without energy recovery, higher extraction efficiency is required. One- 
dimensional computer simulations are used to describe a single pass of the optical wave 
and to optimize the efficiency using a tapered undulator. The tapered undulator represents 
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an advantage when the electron phase acceleration 5 exceeds the deceleration that can be 
obtained without taper in strong optical fields. The single-pass PEL efficiency 7 of a 
tapered undulator is the fraction of the electron beam energy converted to laser light. The 
limit on 7 is determined by the maximum taper rate 8 that can maintain trapped 
electrons over the number of periods N in the undulator [ 6 ]. Gain G is the fractional 
change in power of the optical field in a single pass through the undulator. In weak fields 

without tapering, the PEL amplifier has high gain described by G(r) « (l/9)e^*''^’' 

[12]. The gain is exponential in x along the undulator with growth rate proportional to 
where j = %N{en KVf' pF ! y^mc^ is the dimensionless current density, 
p=2>y.\(f I{A)I ecTxrl is the electron beam density, and F= "area of electron beam'Varea 
of light beam" is the filling factor. 

Table 3.1 shows the proposed parameters of the MW regenerative amplifier 
developed at LANE. 



Parameters 


Proposed Value 


Beam energy E 


100 MeV 


Beam radius 


0.17 mm 


Energy spread Eyly 


0 . 02 % 


Pulse duration 


20 ps 


Pulse frequency 


500 MHz 


Peak current I 


400 A 


Average current / 


0.2 A 


EIndulator parameter K 


1.71 


EIndulator length L 


6 m 


Number of undulator periods N 


300 


EIndulator period 


2.0 cm 


Optical wavelength A 


1.0 pm 



Table 3.1 . Parameters for the MW regenerative amplifier design developed at LANE 
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1, ID Efficiency Optimization 

The proposed MW RAFEL uses a beam of 1 00 Me V electrons with a peak current 
of 400 A, yielding a peak electron beam power of 40 GW. The average current of the 
RAFEL is 0.2 A, so that the average electron beam power is 20 MW. The single-pass 
FEL efficiency r] is the fraction of power extracted from electron beam in one pass 
through the undulator. The RAFEL will feedback on the order of 0.01% to 1% of the 
optical power and require an extraction efficiency of approximately 10% to 15% to 
provide 2 MW to 3 MW of optical power in the infrared (IR). 

In an attempt to achieve the desired efficiency, the initial undulator length is Z = 
4m corresponding Xo N - 200 periods. The dimensionless initial optical field and 

tapered undulator of strength 5 are varied to find the optimum efficiency. The maximum 
extraction efficiency found was 8% which is less than the desired efficiency. 

In an effort to improve the design efficiency, the length of the undulator is 
increased from 4 m to 6 m. The change in undulator length corresponds to dimensionless 
current density of j « 12,000 with N = 300 undulator periods. The values of and d 
are varied to find the optimum efficiency using numerical simulations. The electron 
beam has a Gaussian spread =1.0, and the field strength is =150 with taper rate 

S =260;r starting at time =0.19 along the undulator. Figure 3.2 is the 3D plot of 
efficiency versus and 5 . The maximum efficiency rj^ is 13.3% with gain G = 556. 
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Figure 3.2. 3D plot of Efficiency vs. Optical Field Amplitude 
and Tapered Field Strength. 



25 



The effect of changing away from its optimum value is smaller than changing 5 away 
from its optimum value as shown in Figure 3.3. 




Fig 3.3. Contour plot of Efficiency vs. Optical Field Amplitude 
and Tapered Field Strength. 
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For a better understanding of the physical mechanism underlying the operation of 
this MW RAFEL, examine Figure 3.3 closely in five regions. Region I is the "optimum 
region". As shown in Figure 3.4, the taper turns on at = 0.19 , and approximately half 

of electrons were trapped near resonance = 0 , resulting in an acceptable efficiency of 
7=13.3%. The final gain G^=556 corresponds to feedback of 
/ =1/(G^ + l)=1.8x 10'\ The large optical phase shift (f) in the lower right-hand portion 

of Figure 3.4 represents good optical guiding which is an important part of a high gain, 
high current FEL. 




Figure 3.4. Phase space evolution for the MW RAFEL at optimum efficiency. 
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Region II is the "high tapering rate region". Increasing the tapering rate S beyond 
the optimum value, the efficiency drops to 8.3% due to fewer electrons being trapped near 
resonance = 0 as shown in Figure 3.5. The taper turns on at =0.19 as before. The 

O 

final gain = 347 corresponds to feedback of /=1/(G^ +I)=2.9xl0"\ ^ is 
relatively small which would cause less optical guiding. 



*** FH Phase Spatw Ivolutian'"*** 
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On>l 8>3SOn t)>0.003 




Figure 3.5. Phase space evolution for the MW RAFEL with high tapering rate. 
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Region III is the "low tapering rate region". For 5 below the optimum value, the 
efficiency drops to 6.5% due to insufficient taper rate as shown in Figure 3.6. While many 
electrons are trapped near resonance = 0 , the low taper rate does not extract 
significant energy from them. The taper again turns on at =0.19. The final gain 
Gf = 294 corresponding to feedback of / = 1 l{G^ + 1) = 3 .4 x 1 0’\ 
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Figure 3.6. Phase space evolution for the MW RAFEL with low tapering rate. 
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Region IV is the "high optical field region". Starting with high optical field 
strength , the separatrix does not have a chance to grow and trap electrons, resulting in 
a lower efficiency of 9.5% as shown in Figure 3.7. The taper turns on at = 0.16 . The 
final gain =144 corresponds to feedback of /=l/(Gy + l) = 6.9x 10"^ . 





Figure 3.7. Phase space evolution for the MW RAFEL with high optical field strength. 
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Region V is the "low optical field region". Starting with low optical field strength 
, the initial separatrix is too small so that more electrons leak out from the separatrix, 
resulting in a lower efficiency of 10.9% as shown in Figure 3.8. The taper turns on at 
r^=0.32. The final gain Gf= 25,700 corresponds to feedback of 
+l)=3.9xl0■^ 
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Figure 3.8. Phase space evolution for the MW RAFEL with low optical field strength. 

The above results for efficiency should be considered an over-estimate since the effect of 
diffraction of the optical beam or longitudinal pulse effects are not included. 
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2. 3D Diffraction Optimization 

From Chapter II, the parabolic wave equation can be written as 

a(x,y,z,z)^- j{e"^y (3.1) 

Equation (3.1) governs the dynamics of the optical wave over many optical wavelengths in 
the longitudinal z-dimension, and the operator properly describes the diffraction of the 
optical wave in the transverse (x,_y) directions. A coherent, freely-propagating optical 
wave of radius Wo will spread due to natural diffraction as it traverses the undulator length 
L. The Rayleigh length, nwllX, is the characteristic distance over which the initial 

optical mode area doubles in size. If the optical mode area becomes much larger than the 
electron beam, the optical field will decouple from the electron beam. Therefore, it is 
important that the Rayleigh length is comparable to the undulator length L, and the optical 
mode waist is not too much larger than the electron beam radius rt [13]. 

Three-dimensional simulations are used to describe a single pass of the optical 
wave through the undulator. Initially, the optimal values from one-dimensional simulation 
were used as the inputs for a 3D simulation. The resulting 3D simulation is shown in 
Figure 3.9 with = 150 , d - 260 k, j = 2.6 x lO'*, and the taper turning on at = 0.2 . 

The evolution of the optical mode, |<3f(x,r)| shown in the upper-left window of Figure 3.9, 
is dominated by diffraction after the taper is turned on at = 0.2 . The intense electron 

beam provides guiding over the first one fifth of the undulator, but after tapering begins 
diffraction overcomes this focusing over the rest of the undulator. The top-center window 
shows a cross-sectional view of the optical mode, |a(x,_y)|, at the end of the undulator, 

and this optical mode has a drop in optical power in the center. This drop is due to a loss 
of guiding of the optical mode [14]. The window in the top-right lists the dimensionless 
parameters for this simulation; cr,. and cr_^ correspond to the radial size of the electron 

beam, while (7g and cr^ account for beam quality. 
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FEL WAVEFRONTS 




r, =0.2 77 = 0.0175 J = 2607T Gr=370 

Figure 3.9. 3D diffraction simulation. 



The graph at the right-middle shows the development of the optical phase along the 
undulator length. At =0.2, the optical phase change along the undulator length 

becomes negative and drops rapidly which causes the optical wavefront to be excluded 
from the electron beam decreasing the interaction strength. The lower-right window is the 
plot of the dimensionless optical power growth along the undulator, and the plot of the 
natural logarithm of the single-pass gain along the undulator. At the middle-left shows the 
bunching current in the electron beam development along the undulator, with the end view 
in the center. The bottom-left plot is the electron beam phase velocity evolution along the 
undulator, and at the bottom center shows the electron phase space plot at the end of the 
undulator. Initially the optical guiding confines the optical mode until the taper turns on at 
= 0.2 , then diffraction effects begin to spread the beam. The optical phase begins to 

grow until the taper turns on, then the change in optical phase is negative as the guiding 
effect is lost. The bottom-center window shows that at the end of undulator, no electrons 
are trapped at resonance, and the tapering effect has failed. 
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In order to search for an optimum efficiency using numerical 3D simulations, the 
initial field Go, the taper rate <5, and the time to turn on the taper are varied. The 

optical field amplitude Go is varied from 50 to 275 in steps of 25. For each value of the 
taper field rate 5 is varied from 30;r to 330;r in steps of 2Stu , and for each value of 
and 5 , the turn-on time is varied from 0.1 to 0.5 in steps of 0.02. First we explore the 
MW RAFEL efficiency by varying 5 and at fixed values of Go. 



3D Efficiency Oj3timization (for a.^~ 50) 




Figure 3.10. 3D and Contour plots of Efficiency vs. Tapered Field Strength 
and Taper turn-on time at Gq = 50. 
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Figures 3.10 through 3.14 show the dependence of efficiency// as and 

S are varied. These two parameters can have a significant impact on efficiency 
when moved away from their optimum values. The two graphs of each figure are 
plotted with the same data with the top one is in 3D and the bottom one is in 2D. 



3D Efficiency Optiniiztiiion (for ajj=100) 
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Figure 3.11. 3D and Contour plots of Efficiency vs. Tapered Field Strength 
and Taper tum-on time at Oo= 100. 
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At the initial optical field amplitude of ao = 50, the efficiency peaks at « 6.7% 
at values of ^ = 200;r and = 0.28 . As ao increases, the efficiency peaks at slightly 
higher efficiency «7.0% at lower values of S = \90^ and =0.21 as shown in 
Figure 3.12. 



3D Efficiency Optimization (for a^=175) 
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Figure 3.12. 3D and Contour plots of Efficiency vs. Tapered Field Strength 
and Taper turn-on time at Oo = 175. 
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The peak efficiencies 77 ^^ continue to increase to a optimum efficiency 
'7-2% corresponding to optimum values ao = 225, S = 180;r, and taper turn-on 
time of Tj = 0.18 as shown in Figures 3.13. 



3D Efficiency Opfimization (for 3^=225) 





Contour Plot 




Figure 3.13. 3D and Contour plots of Efficiency vs. Tapered Field Strength 
and Taper turn-on time at Oo = 225. 
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Increasing optical field strength beyond its optimum value, the separatrix does 
not have a chance to grow and trap electrons, resulting in a slightly lower efficiency of 
7.1% as shown in Figure 3.14. 



3D Efllcicncy Opfimization (for a^j=250) 





Contour Plot 




Figure 3.14. 3D and Contour plots of Efficiency vs. Tapered Field Strength 
and Taper tum-on time at Co = 250. 
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Next, explore the effect on efficiency by varying ao and at a fixed value of S . 
As shown in Figure 3.15, starting at a low value of the taper field rate S = S5tz , the 
maximum efficiency obtained was 3.2% with z^ »0.2. This low efficiency is due to 
insufficient taper rate 



3D Efficiency Optimiziifion (for 6=557i:) 
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Contour Plot 




Figure 3.15. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 
and Taper turn-on time aX d = S5n . 
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As S increases, the efficiency increases to « 5.3% at a slightly lower as 
shown in Figure 3.16. The peak efficiency increases to 7^^ ^1.1% corresponding to 
optimum values of (5 = 1 ZOn , Go =225, and an earlier taper turn-on time of = 0. 1 8 . 



3D Efficiency Optimization (for 5=^1057t) 




Figure 3.16. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 
and Taper tum-on time at <5 = 1 05;r . 



40 




As shown in Figure 3.17, when the optimum value of taper rate is used, 5 - 180;^ , 
varying the field Uo has little effect on efficiency, whereas changing the taper turn-on time 
from its optimum value, « 18 , results in a significant efficiency drop. 



3D Efficiency Opf imizution (for S=180 ti:) 




Figure 3.17. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 
and Taper turn-on time at ^ = 180;r . 
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Increasing S beyond its optimum value causes the efficiency to drop significantly 
due to fewer electron being trapped. As shown in Figure 3.18, the efficiency drops to 
only 1.9% at r, = 0.22 with S = ISOtt . 



3D Efficiency Optimization (for 6=280jr) 
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Figure 3.18. 3D and Contour plots of Efficiency vs. Optical Field Amplitude 
and Taper tum-on time at ^ = 2%0n . 
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Now consider the three-dimensional simulation of a single pass of the optical wave 
through the undulator at the optimum values of S = 180;r , = 0.18, and Ug = 225. As 

shown in Figure 3.19, optical guiding occurs and stays with the intense electron beam 
through the undulator as seen in upper-left and middle-left windows. The top-center 
window shows a much tighter cross-sectional end view of the focused optical power than 
seen in Figure 3.9. The bottom-center window shows approximately 50% of the electrons 
are trapped, resulting in optimized gain and efficiency. The middle-right window shows 
the optical phase initially increasing until the taper turns on at « 0.18 , then decreasing 
slowly. The FEL performance in Figure 3.19 is much improved over that of Figure 3.9. In 
steady-state operation, the final gain Gy =616 corresponds to feedback of 

/ = l/(Gy-M) = 1.6xlO-^ 



TEL WAVEFRONTS 




X, =0.18 T1 = 0-072 



I 557*0 * 
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Figure 3.19. 3D diffraction simulation at optimum values 
5 = \ %Qtu , = 0. 1 8 , and ag = 225. 
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IV. CONCLUSIONS 



The development of a speed-of-Iight hard-kill weapon system for military 
applications may represent a significant advancement in technology over present 
conventional kinetic weapon systems. The PEL has the potential for use as a weapon 
system against anti-ship and theatre ballistic missiles. The infinite magazine, rapid 
response, and wavelength tunability at high powers make it suitable and desirable for 
shipboard self-defense and for land-based theater ballistic missile defense. 

Three dimensional computer simulations indicate that the proposed LANL 
Regenerative MW PEL design does not yield the desired extracted efficiency of 77 w 15% . 

The optical field amplitude Qo, taper rate 5 , and taper tum-on time were varied to 
search for the optimum efficiency, resulting in maximum efficiency of only w 7.2% . 
Changing 5 or away from their optimum values has a significant effects on the MW 

RAPEL design efficiency, but changes in Oo do not have much of an effect. The efficiencies 
found here are still an over-estimate since they do not include longitudinal pulse effects. It 
remains to be determined whether a smaller efficiency is acceptable or if 1 5% efficiency is 
actually necessary in the PEL weapon design. If 15% efficiency is necessary, a longer 
undulator, perhaps 8 or 10 meters may be sufficient to reach that goal. 
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